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Abstract. Electron-positron annihilation into hadrons plus an energetic photon from initial state radiation 
allows the hadronic cross-section to be measured over a wide range of energies. The full next-to-leading 
order QED corrections for the cross-section for e'''e~ annihilation into a real tagged photon and a virtual 
photon converting into hadrons are calculated where the tagged photon is radiated off the initial electron or 
positron. This includes virtual and soft photon corrections to the process e'^'e” —>7-1-7* and the emission 
of two real hard photons: e^e“ —> 7 -I- 7 -I- 7*. A Monte Carlo generator has been constructed, which 
incorporates these corrections and simulates the production of two charged pions or muons plus one or 
two photons. Predictions are presented for centre-of-mass energies between 1 and 10 GeV, corresponding 
to the energies of DA^NE, CLEO-C and B-meson factories. 


1 Introduction 

Electroweak precision measurements have become one of 
the central issues in particle physics nowadays. The recent 
measurement of the muon anomalous magnetic moment 
Ufj, = {g — 2)^/2 at BNL [Q led to a new world average, 
differing by 2.6 standard deviations from its theoretical 
Standard Model evaluation [|j. This disagreement, which 
has been taken as an indication of new physics, has trig¬ 
gered a raving and somehow controversial deluge of non¬ 
standard Model explanations. A new measurement with 
an accuracy three times better is under way; this will chal¬ 
lenge even more the theoretical predictions. 

One of the main ingredients in the theoretical predic¬ 
tion for the muon anomalous magnetic moment is the 
hadronic vacuum polarization contribution 0 which is 
moreover responsible for the bulk of the theoretical error. 
It is in turn related via dispersion relations to the cross- 
section for electron-positron annihilation into hadrons, 
o'had = cr(e+e“ ^ hadrons). This quantity plays an im¬ 
portant role also in the evolution of the electromagnetic 
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coupling aQED from the low energy Thompson limit to 
high energies ||^,^. This indeed means that the interpre¬ 
tation of improved measurements at high energy colliders 
such as LEP, the LHC, or the Tevatron depends signifi¬ 
cantly on the precise knowledge of crhad- 

The evaluation of the hadronic vacuum polarization 
contribution to the muon anomalous magnetic moment, 
and even more so to the running of agso, requires the 
measurement of ahad over a wide range of energies. Of 
particular importance for the QED coupling at Mz is the 
low energy region around 2 GeV, where ahad is at present 
poorly determined experimentally and only marginally con¬ 
sistent with the predictions based on pQCD. New efforts 
are therefore mandatory in this direction, which could help 
to either remove or sharpen the discrepancy between the¬ 
oretical prediction and experimental results for (g — 2)^ 
and provide the basis for future more precise high energy 
experiments. 

The feasibility of using tagged photon events at high 
luminosity electron-positron storage rings, such as the (j)- 
factory DA<?NE or H-factories, to measure ahad over a 
wide range of energies has been proposed and studied in 
detail in (see also |^,||]). In this case, the machine 

is operating at a fixed centre-of-mass (cms) energy. Initial 
state radiation (ISR) is used to reduce the effective energy 
and thus the invariant mass of the hadronic system. The 
measurement of the tagged photon energy helps to con- 




2 


G. Rodrigo, H. Czyz, J.H. Kiihn and M. Szopa: Radiative return at NLO 


strain the kinematics, which is of particular importance for 
multimeson final states. In contrast to the conventional en¬ 
ergy scan the systematics of the measurement (e.g. 
normalization, beam energy) have to be taken into ac¬ 
count only once, and not for each individual energy point 
independently. 

Radiation of photons from the hadronic system (final 
state radiation, FSR) should be considered as background 
and can be suppressed by choosing suitable kinematical 
cuts, or controlled by the simulation, once a suitable model 
for this amplitude has been adopted. From studies with 
the generator EVA one finds that selecting events with 
the tagged photons close to the beam axis and well sep¬ 
arated from the hadrons indeed reduces FSR to a rea¬ 
sonable limit. Furthermore, the suppression of FSR over¬ 
comes the problem of its model dependence, which should 
be taken into account in a completely inclusive measure¬ 
ment I p^ . 

Preliminary experimental results using this method 
have been presented recently by the KLOE collaboration 
at DA^NE Large event rates were also ob¬ 

served by the BaBar collaboration |]^ . 

In this paper we consider the full next-to-leading order 
(NLO) QED corrections to ISR in the annihilation process 
e+e“ —> 7 -f hadrons where the photon is observed un¬ 
der a non-vanishing angle relative to the beam direction. 
The virtual and soft photon corrections and the con¬ 
tribution of the emission of a second real hard photon are 
combined to obtain accurate predictions for the exclusive 
channel e“''e“ —> 7r+7r“7 at cms energies of I to 10 GeV, 
corresponding to the energies of DA^NE, CLEO-C and 
R-meson factories. An improved Monte Carlo generator, 
denoted PHOKHARA, includes these terms and will be 
presented in this work. The comparison with the EVA 
Monte Carlo, which simulates the same process at lead¬ 
ing order (LO) and includes additional collinear radiation 
through structure function (SF) techniques, is described. 
Predictions are presented also for the muon pair produc¬ 
tion channel e“'"e“ —> which is also simulated with 

the new generator. 


2 Virtual and soft photon corrections to ISR 

At NLO, the e“'"e“ annihilation process 

e+(pi)-ke"(p 2 ) ^7*(Q)+7(fci) , (1) 

where the virtual photon converts into a hadronic final 
state, 7*((5) ^ hadrons, and the real one is emitted from 
the initial state, receives contributions from one-loop cor¬ 
rections (see Fig. P) and from the emission of a second real 
photon (see Fig. g). 

After renormalization the one-loop matrix elements 
still contain infrared divergences. These are cancelled by 
adding the contribution where a second photon has been 
emitted from the initial state. This rate is integrated ana¬ 
lytically in phase space up to an energy cutoff < Wyfs 
far below The sum |l6| is finite; however, it depends 
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Fig. 1. One-loop corrections to initial state radiation in 
e'*'e“ ^7-1- hadrons. 


on this soft photon cutoff. The contribution from the emis¬ 
sion of a second photon with energy Ej > w^/s completes 
the calculation and cancels this dependence. 

In order to facilitate the extension of the Monte Carlo 
simulation to different hadronic exclusive channels the dif¬ 
ferential rate is cast into the product of a leptonic and a 
hadronic tensor and the corresponding factorized phase 
space: 


da = -^Li,,,Hf^^d<p2{Pi,P2]Q,ki)d^n{Q;qi,-,qn)‘^^ , 
Zs Ztt 

( 2 ) 

where d<Pn{Q', qi, ■, q-a) denotes the hadronic n-body phase 
space including all statistical factors and is the invari¬ 
ant mass of the hadronic system. 

The physics of the hadronic system, whose description 
is model-dependent, enters only through the hadronic ten¬ 
sor 


^ j^+ ^ (3) 

where the hadronic current has to be parametrized through 
form factors. For two charged pions in the final state, the 
current 


JL = {q^+ - q^-Y (4) 

is determined by only one function, the pion form factor 
F 2 Tr [P^ . The hadronic current for four pions exhibits a 
more complicated structure and has been discussed in [Q . 

The leptonic tensor, which describes the NLO virtual 
and soft QED corrections to initial state radiation in e+e“ 
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annihilation, has the following general form: 
( 47 ra)^ 
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where yi = 2 ki ■ pi/s. The scalar coefficients Oy and o_i 
allow the following expansion 
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where o|°^ give the LO contribution. The NLO coefficients 

o^j^ and were calculated in |Q (see also [^) for the 
case where the tagged photon is far from the collinear 
region. 

The expected soft and collinear behaviour of the lep- 
tonic tensor is manifest when the following expression is 
used: 
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where stands for the LO leptonic tensor. The first 
term, log( 4 t(;^)[l + log(TO^)] where = m^/s, contains 
the dependence on the soft photon cutoff lu, which has to 
be cancelled against the contribution from hard radiation. 
The next three terms, also proportional to the LO leptonic 
tensor, represent the QED corrections in leading log ap¬ 
proximation with the typical logarithmic dependence on 
the electron mass. The tensor finally contains the 
subleading QED corrections. 

The factor 1/(1 — Syp) accounts for the vacuum polar¬ 
ization corrections in the virtual photon line. This multi¬ 
plicative correction can be approximately reabsorbed by 
a proper choice of the running coupling constant. In the 
present version of the Monte Carlo generator one can 
choose to include or not the contribution from the real 
part of the lowest order leptonic loops p^: 
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Fig. 2. Emission of two real photons from the initial state in 
e^e“ annihilation into hadrons. 


3 Emission of two hard photons 

In this section, the calculation of the matrix elements for 
the emission of two real photons from the initial state 


log(4r(;^)[l -I- log(m^)] 


e+(pi)-he (P 2 ) ^ 7*(Q)+ 7 (^ 1 )+7(fe) , 


( 11 ) 


is presented (see Eig. ||). 

We follow the helicity arnplitude method with the con¬ 
ventions introduced by |^,^. The Weyl representation 
for fermions is used where the Dirac matrices 




0 cr(/ 

7 ? 0 


M = 0,1,2,3 


( 12 ) 


are given in terms of the unit 2x2 matrix / and the 
Pauli matrices cri,z = 1 , 2 , 3 , with = (/, ±cri). The 
contraction of any four-vector with the 7^^ matrices 
has the form 


0 a+ 
a~ 0 


4 = 07*7'" = 

where the 2x2 matrices are given by 


= 


a? =F 0? =F(a^ - io?) 

7 ^{a^+io?) a°±a^ 


(13) 


( 14 ) 


The helicity spinors u and v for a particle and an 
antiparticle of four-momentum p — (E, p) and helicity 
A /2 = ± 1/2 are given by 


V)- ( x(p,A)\ _ 

- f-W^+M pI X(P: -?) 

A/E^^x(p,-A) 
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( 15 ) 


The helicity eigenstates x(P, "^) can be expressed in terms 
of the polar and azimuthal angles of the momentum vector 
p as 


x(p,+i) = 
x(p,-i) = 


cos(0/2)A 
sin (0/2 )) 


-e-*'^sin 


(^/2)\ 
cos(0/2 )J 


This routine can be easily replaced by a user if necessary. 


(16) 
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Finally, complex polarization vectors in the helicity basis 
are defined for the real photons: 

= =p) = —^(0,±cosdiCOS(/)i +*sin(/)i, 

v2 

± COS Oi sin (pi — i cos pi , =F sin 9 i) , 

( 17 ) 


with i = 1 , 2 . 

The complete amplitude can be written in the follow¬ 
ing form 

M (Ae+, Ae- , Ai, A2) = v\{\f.+ )A (Ai, A2) Ul{Xe- ) 

+ vlj{Xe+)B {\i, X 2 ) UII (Xf.-) , {18) 

where A (Ai, A2) and B (Ai, A2) are 2x2 matrices defined 
in the appendix with the matrix elements (Ai, A2) and 
Bij{Xi, X2). It simplifies even further, when calculated in 
the electron-positron cms frame (the z-axis was chosen 
along the positron momentum): 


M (-I-, -1-, Ai, A2) = TOe {A 2,2 (Ai, A2) — -82,2 (Ai, A2)) 
M (—, —, Ai, A2) = nie {Ai^i (Ai, A2) — Bi i (Ai, A2)) 
M (-I-, —, Ai, A2) = — {E + |p|) ^2,1 (Ai, A2) 


+ - 


mz 


782,1 (Ai, A2) , 


E + ipl 

M {—, -|-, Ai, A2) = (8 -f |p|) 81,2 (Ai, A2) 

7^1,2 (Ai, A2) 


8 - 


( 19 ) 


From the explicit form of the matrices A(Ai,A2) and 
8 (Ai, A2), it is clear that some factors appear repeatedly 
in different amplitudes. In order to speed up the numerical 
computation, the amplitudes are decomposed into these 
factors, which are used as building blocks for all the di¬ 
agrams. Then, a polarized matrix element is calculated 
numerically for a given set of external particle momenta 
in a fixed reference frame, e.g. the cms of the initial parti¬ 
cles where the initial momenta are parallel to the z axis. 
The result is squared and the sum over polarizations is 
performed. 

As a test, the square of the matrix element averaged 
over initial particle polarization has also been calculated 
using the standard trace technique and tested numerically 
against the helicity method result. Perfect agreement be¬ 
tween the two approaches is found. Both matrix elements 
vanish if the photon polarization vectors are replaced by 
their four-momenta, and thus are tested for gauge invari¬ 
ance. 


the previous versions of the EVA 0 Monte Carlo gen¬ 
erator. The program exhibits a modular structure, which 
preserves the factorization of the initial state QED correc¬ 
tions. The simulation of other exclusive hadronic channels 
can therefore be easily included with the simple replace¬ 
ment of the current(s) of the existing modes, and the cor¬ 
responding multiparticle hadronic phase space. The sim¬ 
ulation of the four-pion channel H will be incorporated 
soon, as well as other multihadron final states. 

The program provides predictions either at LO or at 
NLO. In the former case only single photon events are 
generated. In the latter, both events with one or two pho¬ 
tons are generated at random. For simplicity, FSR is not 
considered in the new generator, which however can be 
estimated from EVA . 

Single photon events are generated following the same 
procedure as in EVA. The generation of two-photon events 
proceeds as follows. First, polar and azimuthal angles of 
the two photons are generated. One of the polar angles is 
generated within the given angular cuts, the other is gener¬ 
ated unbounded. In this way, the photon angular cuts are 
automatically fulfilled and a higher generation efficiency 
is achieved. Both photons are nevertheless symmetrically 
generated. Then, the hadronic invariant mass is gen¬ 
erated following the resonant distribution of the hadronic 
current, its maximum being determined by eq.(p^ in Ap¬ 
pendix where one of the photon energies is set to the 
soft photon cutoff and the other is set to the minimal de¬ 
tection photon energy 8™®”. Next, the photon energies are 
generated. If only one of the photons passes the angular 
cuts, its energy is forced to be larger than 8("™. The other 
photon energy is calculated according to eq. (^ . Other¬ 
wise, if both photons pass the angular cuts the minimal en¬ 
ergy of one of them is fixed to w, with equal probability for 
both photons, and the other is calculated through eq.(^. 
Finally, the hadron momenta are generated in the rest 
frame and then boosted to the e^e“ cms. Further angu¬ 
lar cuts or other kinds of constraints are imposed after 
generation. For more details, see Appendix H. 


5 Tests and results 

5.1 Gauss integration versus MC 

To test the technical precision of the single photon gener¬ 
ation, a FORTRAN program was written, which performs 
the two-dimensional integral that remains after the inte¬ 
gration over the pion angles, and the photon azimuthal 
angle has been performed analytically with the help of 
the relation 


4 Monte Carlo simulation 

On the basis of these results a Monte Carlo generator, 
called PHOKHARA, has been built to simulate the pro¬ 
duction of two charged pions together with one or two hard 
photons; it includes virtual and soft photon corrections 
to the emission of one single real photon. It supersedes 


I Jf(j^)* d<p2iQ;qi,q2) = 

^ (Qi.Q. - q^rQ") B(Q^) . ( 20 ) 

For the remaining numerical integration, the integra¬ 
tion region was sliced into an appropriate number of subin¬ 
tervals (typically 100 to 200 ), and 8-point Gauss quadra- 
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Fig. 3. Relative difference between the single photon cross 
section obtained from the Monte Carlo generator and 

Gauss quadrature (see text for details). 


Fig. 4. Comparison of the differential distribution for dif¬ 
ferent values of the soft photon cutoff: w = 10“®, 10“'* and 
10"®, at y/s = 1.02 GeV. 


ture was used in each of the subintervals. This leads to a 
relative accuracy of 10"^°. 

The test was performed for a photon polar angle B-y 
between 10° and 170° and for 4m^ < < 2Ecm{Ecm — 

Ernin) ^ ^ GeV and = 10 MeV. To 

study contributions from different ranges of 9^ and 
the above-mentioned intervals of and cos(0.y) were di¬ 
vided into ten equally spaced parts. The integrals were 
performed first within the whole range of 9^, and ten 
subintervals in separately and subsequently within the 
whole range of and ten subintervals of cos(0.y). From 
Fig.0 it is clear that a technical precision of the order of 
10“^ was achieved. The error bars indicate one standard 
deviation of the Monte Carlo generator, which performs 
the five-dimensional integral. 


5.2 Soft photon cutoff independence 

The full NLO calculation consists of two complementary 
contributions, the virtual and soft corrections presented in 
section ^ and the hard correction described in section ||. 
The former depends logarithmically on the soft photon 
cutoff w, see eq.(^. The second, after numerical integra¬ 
tion in phase space exhibits the same behaviour, so that 
their sum must be independent of w. 

However, a particular value of w has to be chosen for 
the generation. To be valid, the soft photon approximation 
requires w to be small. On the other hand a very small 
value of w could even produce unphysical negative weights 
for the generated events. The particular value of w chosen 
to perform the Monte Carlo generation should therefore 
arise from a compromise between these two conditions. 
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Fig. 5. Comparison of the differential distribution for dif¬ 
ferent values of the soft photon cutoff: w = 10~®, 10"'* and 
10"®, at = 10.6 GeV. 
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Table 1. Kinematical cuts applied at different cms energies: 
minimal energy of the tagged photon (-E 7 ), angular cuts on the 
tagged photon {9^) and the pions {6^), and minimal invariant 
mass of the hadrons plus the tagged photon 

yi=1.02GeV 4GeV 10.6 GeV 
(GeV) 0.01 0.1 1 

Oj (degrees) [5,21] [10,170] [25,155] 

9^ (degrees) [55,125] [20,160] [30,150] 

(GeV")0^1^90 


Table 3. Total cross section (nb) for the process e''"e“ —> 
7 r"'" 7 r ”7 at LO, NLO(l) and in the collinear approximation via 
structure functions (SF) with the cuts from Table Only ini¬ 
tial state radiation. NLO(2) gives the NLO result with the same 
cuts as NLO(l) but for the hadron-photon invariant mass un¬ 
bounded. 

=1.02 GeV 4 GeV 10.6 GeV 
Born 2.1361 (4) 0.12979 (3) 0.011350 (3) 

SF 2.0192 (4) 0.12439 (5) 0.010526 (3) 

NLO (1) 2.0332 (5) 0.12526 (5) 0.010565 (4) 

NLO (2) 2.4126 (7) 0.14891 (9) 0.012158 (9) 


Table 2. Total cross section (nb) for the process e'’"e“ —> 
7 r^ 7 r “7 at NLO for different values of the soft photon cutoff. 
Only initial state radiation. Cuts from Table ^ 


w 

ys = 1.02 GeV 

4 GeV 

10.6 GeV 

10 “^ 

2.0324 (4) 

0.12524 (5) 

0.010564 (4) 

10 “^ 

2.0332 (5) 

0.12526 (5) 

0.010565 (4) 

10 “® 

2.0333 (5) 

0.12527 (5) 

0.010565 (5) 


In this section we show that the result from the gen¬ 
erator is indeed independent from the soft photon cutoff 
tc, within the error of the numerical integration. Further¬ 
more, we try to determine the value of w that optimizes 
the event generation. 

The tests have been performed for different cms en¬ 
ergies, from 1 to 10 GeV, corresponding to the energies 
of DAl^NE, CLEO-C and i?-meson factories. Kinematical 
cuts have been applied as listed in Table |] and will be used 
for the rest of the paper. They are related to those of the 
experiments for which we present our predictions and at 
the same time allow final state radiation to be controlled, 
one of the key points of the radiative return method. A 
minimal energy E™*” is required for the tagged photon. 
Different cuts are chosen for the polar angle of the tagged 
photon 9^ at different energies. At low energies, the pi¬ 
ons are constrained to be well separated from the pho¬ 
tons to suppress final state radiation. At high energies, 
the observed photon and the pions are mainly produced 
back to back; the suppression of the final state radiation 
is therefore naturally accomplished. Furthermore, a mini¬ 
mal invariant mass of the hadronic system plus the tagged 
photon, is required. The reason for this last kine¬ 

matical cut will be discussed later. When events with two 
photons are simulated we require at least one of the pho¬ 
tons to pass the cuts. 

Table || presents the total cross section calculated for 
several values of the soft photon cutoff at three different 
cms energies for the kinematical cuts from Table ^ The 
excellent agreement confirms the ic-independence of the 
result. 

Figures ^ and || show the rc-independence of the 
distribution at ^/s = 1.02 GeV and 10.6 GeV cms ener¬ 
gies respectively. Even if the analysis of the total cross 
section (Table g) suggests that the choice of ic = 10“^ is 
as good as w = 10“^, small differences in the differential 
cross section are found for high values of Q^. Similar com¬ 
parisons between the differential cross section calculated 


Table 4. Total cross section (nb) for the process e"'"e“ —» 
7 r"'" 7 r ”7 at =1.02 GeV in NLO and in the collinear approx¬ 
imation (SF) as a function of the cut on the invariant mass of 
the hadron -|- tagged photon Only initial state radi¬ 

ation. Minimal energy of the tagged photon and angular cuts 
from Table 


(GeV^) 

SF 

NLO 

0.1 

2.4127(18) 

2.4132(8) 

0.2 

2.4126(18) 

2.4131(8) 

0.3 

2.4124(18) 

2.4127(8) 

0.4 

2.4098(18) 

2.4096(8) 

0.5 

2.3949(18) 

2.3953(8) 

0.6 

2.3425(16) 

2.3455(8) 

0.7 

2.2449(11) 

2.2543(8) 

0.8 

2.1387(9) 

2.1533(8) 

0.9 

2.0198(8) 

2.0334(8) 

0.95 

1.9437(8) 

1.9522(8) 

0.99 

1.8573(8) 

1.8559(8) 


for w = 10“"^ and w = 10“® show a perfect agreement 
within the statistical errors, well below the 0.1% level. As 
a result we chose w = 10“^ for the soft photon cutoff. 


5.3 Comparison of NLO results with the structure 
function approximation and an estimation of the 
theoretical error 


The original and default version of EVA g, simulating 
the process e+e” ^ 7r+7r“7 at LO, allowed for additional 
initial state radiation of soft an d collinear photons by 

I. By convoluting the 


the structure function method | ^ 


Born cross section with a given SE distribution, the soft 
photons are resummed to all orders in perturbation theory 
and large logarithms of collinear origin, L = log{s/ml), 
are taken into account up to two-loop approximation. The 
NLO result, being a complete one-loop result, contains 
these logarithms in order a and the additional subleading 
terms, which of course are not taken into account within 
the SF method. The subleading terms from virtual and 
soft corrections were calculated in [|^ and are included in 
the present NLO generator. 

In the SF approach, the additional emission of collinear 
photons reduces the effective cms energy of the collision. 
Ini, a minimal invariant mass of the observed particles. 
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Fig. 6. Differential cross section for the process e''"e —*■ 

7r'^7r“7 at NLO for ^/s = 1.02 GeV. Only initial state ra¬ 
diation. The relative size of the correction to the LO result 
(dashed) is shown in the small inset. The cuts are: 5° < 
9^ < 21°, 55° < Ott < 125°, the energy of the tagged pho¬ 
ton Ej > 0.01 GeV and the invariant mass of the detected 
particles in the final state > 0.9 GeV^ for NLO(l) 

(solid). NLO(2) (dotted) obtained without the last cut. 


Fig. 8 . Differential cross section for the process e''"e“ —> 
7r"'"7r~7 at NLO for ^/s = 10.6 GeV. Only initial state ra¬ 
diation. The relative size of the correction to the LO result 
(dashed) is shown in the small inset. The cuts are: 25° < 
9j < 155°, 30° < Btt < 150°, the energy of the tagged photon 
E~, > 1 GeV and the invariant mass of the detected particles in 
the final state > 90 GeV^ for NLO(l) (solid). NLO(2) 

(dotted) obtained without the last cut. 
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Fig. 7. Gomparison between the collinear approximation by 
structure functions and the fixed order NLO result. Cuts from 
Table ^ for the lower figure. In the upper figure, same cuts as 
below, with the addition of the symmetric photon configuration 
159° <ej < 175°. 
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Fig. 9. Comparison between the collinear approximation by 
structure functions and the fixed order NLO result. Cuts from 


Table 


hadrons plus tagged photon, was required in order to re¬ 
duce the kinematic distortion of the events. To perform 
the comparison between EVA and the present program 
a similar cut is introduced in the NLO calculation. For 
one-photon events the invariant mass of the hadrons and 
the emitted photon is equal to the total cms energy and 
the requirement is trivially fulfilled. However, when two- 
photon events are generated, energy and production an¬ 
gles of both photons are compared with the cuts listed in 
Table |^. If one of the photons fulfils both requirements, 
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Fig. 10. Differential cross section for the process e'’'e ^ 

7r^7r“7 at NLO for ^/s = 4 GeV. Only initial state radiation. 
The relative size of the correction to the LO result (dashed) is 
shown in the small inset. The cuts are: 10° < 9-y < 170°, 20° < 
07 r < 160°, the energy of the tagged photon E.y > 0.1 GeV 
and the invariant mass of the detected particles in the final 
state > 12 GeV^ for NLO(l) (solid). NLO(2) (dotted) 

obtained without the last cut. 



e" (GeV^) 


Fig. 11. Comparison between the collinear approximation by 
structure functions and the fixed order NLO result. Cuts from 
Table |^. 


its common invariant mass with the hadrons is calculated. 
In other words, we require at least one of the photons to 
pass all the cuts, including the one on its common invari¬ 
ant mass with the hadrons. The probability of both pho¬ 
tons passing all the cuts becomes very small when this 
invariant mass cut is close to the total cms energy. 

Table || gives the total cross section calculated at LO 
and NLO for the previous kinematical cuts. The soft pho¬ 
ton cutoff is fixed at w = 10“^. For comparison, the total 
cross section derived from EVA with emission of collinear 
photons by structure functions [^,|2| is presented. Two 


NLO predictions are shown. The first one, NLO(l), which 
can be compared with the SF result, includes the cut on 
the invariant mass of the hadrons plus photon. The second 
one, NLO(2), is obtained without this cut. The results of 
EVA and those denoted NLO(l) for the total cross section 
are in agreement to within 0.7%. Both of them are clearly 
sensitive to the cut on This cut dependence is 

displayed in Table Remarkably enough, the typical dif¬ 
ference between the results of the two programs is clearly 
less than 0.5% for most of the entries. 

Figure ^ presents our NLO predictions for the differ¬ 
ential cross section of the process e+e" ^ 7r+7r“7 as a 
function of , the invariant mass of the hadronic system, 
at DA^NE energies, ^/s = 1.02 GeV, with the same kine¬ 
matical conditions as before. For comparison, also the LO 
prediction (without collinear emission) is shown, as well as 
the relative correction. Notice that although the NLO(l) 
prediction for the total cross section differs from the LO 
result by roughly 5% (see Table |3|), the distribution 
shows corrections up to ±15% at very low or high values 
of 

As already shown in Table ^ for the total cross sec¬ 
tion, the invariant mass of the hadron ± tagged photon 
system of a sizeable fraction of events lies below the cut 
of 0.9 GeV^. These are events where a second hard pho¬ 
ton must be present. Including these events leads to the 
distribution denoted NLO(2) in Fig. ||. 

In contrast to the typically 5% difference between LO 
and NLO(l) predictions, only a 0.6% difference is found 
between the NLO(l) prediction and the SF collinear ap¬ 
proximation as implemented in EVA, this difference being 
higher at small and lower at high see Fig. As one 
can see from Fig. ^ the size and sign of the NLO correc¬ 
tions do depend on the particular choice of the experimen¬ 
tal cuts. Hence only using a Monte Carlo event generator 
one can realistically compare theoretical predictions with 
experiment and extract R{s) from the data. The difference 
between Figs. ^ and 0b arises from the (small) subset of 
events with two photons, which both fulfil the angular 
and energy cuts and thus enter only once in the sample of 
Fig. 0b. 

Results at 10.6 GeV cms energy are presented in Figs. | 
and 0. In this case a NLO(l) correction to the LO result 
of 7.5% is found, almost independent of Q^, the differ¬ 
ence between the NLO(l) prediction and the SF collinear 
approximation being smaller than 0.5% and also almost 
independent of Q^. 

Finally, predictions for y/s = 4 GeV are presented in 
Figs. 1^ and 0 

To estimate the systematic uncertainty of the program, 
we observe that leading logarithmic two-loop 0{a^) cor¬ 
rections and the associate real emission are not included. 
The difference between the LO and the NLO(l) results 
was expected to be of the order of |(Q:/7r) log (s/rUg) « 
5% at y/s = 1 GeV. Indeed one observes (see Fig. 0) 
values typically of this magnitude with maximal devia¬ 
tions close to 10%. From naive exponentiation one expects 
5 (|(a/7r) log(s/TOg))^ « 0.1-0.2% for the leading pho- 












































G. Rodrigo, H. Czyz, J.H. Kiihn and M. Szopa: Radiative return at NLO 


9 


tonic next-to-next-to-leading order (NNLO) terms, which 
are ignored in the present program. 

Another type of o? corrections originates from fermion 
(dominantly electron) loop insertions in the one-loop vir¬ 
tual corrections considered in 116 . These are conveniently 
combined with the emission of (mainly soft and collinear) 
real fermions, i.e. with the process e“'"e“ ^ 77*//, where 
the collinear // pair is mostly within the beam pipe and 
thus undetected. Individually these are corrections of or¬ 
der (a/7r)^(log(s/TOg))^. When combined, the log^(s/mg) 
terms cancel and the remaining i(Q;/7r)^(log(s/TOg))^ term 
can be largely absorbed by choosing a running coupling 
a(l GeV) in the NLO terms and thus amount to less than 
0.1%. Adding these contributions linearly, one estimates 
a 0.3% uncertainty [|. 

Multiphoton emission is not included in the present 
program. In an inclusive treatment and for tight cuts on 
the photon energy, this can in principle be included through 
exponentiation. For the cuts on of 0.9 GeV^ pro¬ 

posed originally in |^, and adopted through most of this 
paper, the difference between the exponentiated form and 
the fixed order treatment (see e.g. eqs. 17 and 19 of |p^ ) 
amounts to roughly 0.7%, with smaller values for less re¬ 
strictive cuts. For a cut of , _ at 0.8 GeV^ or even 

TT^TT 7 

0.7 GeV^, which seems preferable from these considera¬ 
tions, we expect a difference of 0.4% and even below 0.3% 
for the second choice. In total we therefore estimate a sys¬ 
tematic uncertainty from ISR of around 0.5% in the total 


cross section, once loose cuts on M/ 


are adopted. 


FSR can be controlled by suitable cuts or corrected with 
the Monte Garlo program EVA. 


5.4 Muon pair production 


Inclusion of muon production in the program is straight¬ 
forward. The results for the total cross section are listed 
in Table P, the differential distributions for the three cms 
energies in Figs. The radiative muon cross sec¬ 

tion could be used for a calibration of the pion yield. 
A number of radiative corrections are expected to can¬ 
cel in the ratio. For this reason we consider the ratio be¬ 
tween the pion and the muon yields, after dividing the 
former by |i7((5^)P(I — the latter by 4(1 -|- 

2 m^/(5^)Y^l — 4 m 2 /Q 2 _ pjg^ p^ a we consider the full 

angular range for pions and muons, with between 5° 
and 21°. Clearly all radiative corrections and kinematic 
effects disappear, up to statistical fluctuations, in the lead¬ 
ing order program as well as after inclusion of the NLO 
corrections. 

In Fig. ^0 an additional cut on pion and muon an¬ 
gles has been imposed. As demonstrated in Fig. p^, the 
ratio differs from unity once (identical) angular cuts are 
imposed on pions and muons, a consequence of their dif¬ 
ferent respective angular distribution. To derive the pion 


^ Note that the reaction e''"e“ —> 7 *// leads to significantly 
larger effects However, this process does not contribute 

to events with a tagged photon. 


Table 5. Total cross section (nb) for initial state radiation in 
the process e^e“ —> at LO, NLO (1) and NLO (2) with 

the cuts from Table [iJ, the pions being replaced by muons. 

=1.02 GeV 4 GeV 10.6 GeV 
Born 0.8243(5) 0.4690(6) 0.003088(6) 

NLO (1) 0.7587(5) 0.4449(6) 0.002865(6) 

NLO (2) 0.8338(7) 0.4874(14) 0.00321(6) 




I 

+ 

t 


D 

■« 



o' (GeV^) 


Fig. 12. Differential cross section for the process e''"e —> 

at NLO for = 1.02 GeV. Only initial state radiation. 
The relative size of the correction to the LO result (dashed) 
is shown in the small inset. Same cuts as in figure ^ with the 
pions replaced by muons. 


form factor from the ratio between pion and muon yields, 
this effect has to be incorporated. However, the correction 
function shown in Fig. & is independent from the form 
factor, and hence universal. It can be obtained from the 
present program in a model-independent way (ignoring 
FSR for the moment). 


6 Conclusions 

The radiative return with tagged photons offers a unique 
opportunity for a measurement of cr(e+e“ —> hadrons) 
over a wide range of energies. The reduction of the cross 
section by the additional factor a/ir is easily compensated 
by the high luminosity of the new e“'"e“ colliders, specifi¬ 
cally the /)-, charm- and B-factories at Frascati, Cornell, 
Stanford and KEK. 

In the present work we presented a Monte Carlo gen¬ 
erator that simulates this reaction to next-to-leading ac¬ 
curacy. The current version includes initial state radia¬ 
tion only and is limited to 7r'''7r“7(7) and as fi¬ 

nal states. The uncertainty from unaccounted higher order 
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Fig. 13. Differential cross section for the process e'’'e“ ^ 
fi~'y at NLO for ^/s = 4 GeV. Only initial state radiation. 
The relative size of the correction to the LO result (dashed) 
is shown in the small inset. Same cuts as in Fig. with the 
pious replaced by muons. 



Fig. 14. Differential cross section for the process e'’"e ^ 

at NLO for = 10.6 GeV. Only initial state radiation. 
The relative size of the correction to the LO result (dashed) 
is shown in the small inset. Same cuts as in Fig. ^ with the 
pions replaced by muons. 




Fig. 15. Ratio between pion and muon yields, after dividing 
through their respective R-ratio. Fig. a: no cuts on pion and 
muon angles. Fig. b: with angular cuts on pion and muon an¬ 
gles. 

ISR is estimated at around 0.5%. The dominant FSR con¬ 
tribution can be deduced from the earlier program EVA. 
Additional hadronic modes can be easily implemented in 
the present program. The modes with three and four pions 
are in preparation. 
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A Helicity amplitudes 


In this appendix the full set of helicity amplitudes for the 
diagrams of Fig. || is given. Notation and calculation pro¬ 
cedure were outlined in section which follows Refs. 

As an example, consider the amplitude of the first 
Feynman diagram in Fig. 

Ml = = ^v{pi)f {ki)[}^i -^1 -I-me] 

X f{k2)[Hi + ^2 -#i + rn^]Yu[p2)Jp. 

X _i_ . (21) 

(2fci • pi){2ki ■ k2 — 2ki ■ Pi — 2k2 ■ pi) 

Using the Dirac equation, after some algebra, the following 
expression is obtained 


where 

Ni = ki ■ Pi + k 2 ■ Pi - ki ■ k 2 , 

N2 = ki ■ p2 + k2 ■ P2 - ki ■ k2 , 
ai = £*{k 2 )~k 2 - 2 e*{k 2 ) ■ pi , 

02 = e*{k 2 )'^k 2 — 2 e*{k 2 ) ■ pi , 
a^ = 2 p 2 - J - Q^J~ , 

04 = 2 p 2 ■ J — Q~J'^ , 

05 = e*{ki)~k^ - 2 £*{ki) ■ pi , 

06 = e*(fci)+fcC - 2£:*(fci) • Pi , 

07 = 2 £*{k 2 ) ■ P2 - , 

08 = 2e*(fc2) ■P2 - fc^e*(fc2)+ , 

09 = 2e*(fci) • p2 - k'^e*{ki)~ , 
oio = 2£*(fci) ■p2 - fc]“e*(fci)+ , 
dll — J Q ‘^Pl * J ; 

012 = J~^Q~ - 2pi • J . 

( 24 ) 

The current is defined by the eq.(U) for the tt+tt” 
final state, while for the one it is defined as follows: 

A2) = ieu(( 7 i, Ai) 7 '"w(( 72 , A2) , ( 25 ) 

where <71, Ai (92, A2) are four-momentum and helicity of 

p- {p+). 

B Phase space 


^ (ki) kt - 2£*{ki) -pi] 

X e*(fc 2 )“[ 2 p 2 • J - Q~^J~]ui{p 2 ) 

+ vlj{pi)[£*{ki)^ki -2£*{ki) -pi] 

X e*(fc 2 )^[ 2 p 2 • J - Q~ J~^]uii{p 2 )^ 

X _ - _ . (22) 

(2fci • pi)(2fci • k2 — 2fci • Pi — 2k2 ■ Pi) 

Similar expressions are found for the amplitudes of the 
other diagrams and the complete matrix element can be 
written in a simple form (|l8|), where the matrices A (Ai, A2) 
and B (Ai, A2) are defined as: 


A = - 


B = - 


e" 

W 


0103 


(k2 ■pi){ki ■P 2 ) 
3 


e 

W 


0204 


0503 


(k2 ■pi)Ni 

0109 


0507 


(fci • pi)A^i 
01107 ^ 

{k2 ■P2)N2 

0504 


(fci • Pi){k2 ■ P2) 
01109 


{ki 


■P2)N2 

0608 


{k2-pi)Ni {ki-pi)Ni {ki ■ pi){k2 ■ P2) 


The generation of the multiparticle phase space is based 
on the following Lorentz-invariant representation: 

d<Prn+niPl,P2; ki, ■, km, qi, ■, qn) = 

dQ'^ 

dd>m{pi,P2;Q,ki,-,km)dd>niQ;qi,-,qn)-z - , (26) 

ZTT 

where pi and p 2 are the four-momenta of the initial par¬ 
ticles, ki... km are the four momenta of the emitted pho¬ 
tons and qi... qn, with Q = ^ qi, label the four-momenta 
of the hadrons. 

When two pions are produced in the final state, the 
hadronic part of phase space is given by 


d<p2{Q;qi,q2) 



327r2 


df 2 , 


(27) 


where df 2 is the solid angle of one of the pions at, for 
instance, the rest frame. 

One single photon emission is described by the corre¬ 
sponding leptonic part of phase space 


d‘d>2{pi,P2',Q,ki) 


1-g^ 

327r2 




(28) 


02010 


(fc2 •Pl)(fcl ■P2) 


O12O8 

{k2 ■P2)N2 


012010 A 
(fci-p2)A^2/ 

( 23 ) 


with q^ = Q^/s and df2i is the solid angle of the emitted 
photon at the e+e“ rest frame. The polar angle 6*1 is de¬ 
fined with respect to the positron momentum pi. In order 
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to make the Monte Carlo generation more efficient, the 
following substitution is performed: 

„ 1 , / ^ N 1 . l + / 3 cos 0 i 

cos6»i = - tanh (/3 h) , = 7^ log -- , ( 29 ) 

p Zfj I — fj cos 0 i 

with /3 = — Amljs, which accounts for the collinear 

emission peaks 


The collinear peaks are then flattened with the help of 
eq. (1^ ) , with one change of variables for each photon polar 
angle. The remaining soft peak, wi —> w, is reabsorbed 
with the following substitution 


Wi 


1 -g^ 

2 + e-“i ’ 


ui = log 


Wl 

1 — — 2wi 


( 37 ) 


d cos 01 
1 — cos^ 01 


( 30 ) 


Then, the azimuthal angle and the new variable ti are 
generated flat. 

Consider now the emission of two real photons. In 
the cms of the initial particles, the four-momenta of the 
positron, the electron and the two emitted photons are 
given by 


or 


dwi dui 

wi(l-g2_2wi) “ 1 - ’ 

where the new variable ui is generated flat. 
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